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ABSTRACT 

The application of Mitrovic's Method enables the control systems 
engineer to analyze and synthesize control systems as a function of two 
parameters. The feedback compensation of sampled-data control systems 
is thoroughly investigated as a multi-parameter problem, with emphasis 
placed on the design of the compensation in the Mitrovic parameter plane. 
The minimization of the control system settling time is shown to be parti- 
cularly adapted to the parameter plane. Several examples are included 
which demonstrates this usefulness of the parameter plane in solving the 
minimum settling time problem. A stability relationship is devised for 
digital computer determination of dynamic system stability which permits 
the formulation of an n-dimensioned parameter space. Applications of this 
digital creation of an n-dimensioned space are considered. This ability 
to design feedback compensation for the sampled-data system as a function 
of more than one parameter is shown to lead to greater accuracy in closed 
loop root positioning; and in certain cases, to permit the design con- 


sideration of severa! constraints simultaneously. 
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CHAPTER J 
GENERA], CONCEPTS AND THEORY 
1.1 Introduction. 

The design of multivariable feedback compensation for a feedback 
control system, linear or sampled-data, has heretofore been accomplished 
largely by one of two methods. Either repeated calculations using the 
classical one parameter design techniques, such as root locus, has been 
used; or, if a more accurate study is required, the designer could resort 
to analog simulation and a “trial and error" approach to the design pro- 
blem. Dusan Mitrovic, {3}, in his early work with the algebraic methods 
for locating the roots of a polynomial introduced a new analysis technique 
that permitted analysis of the root location in terms of two variable co- 
efficients. Later, Siljak, [4] and [5], devised additional algebraic 
methods whereby two parameters could be selected from a portion of any, 
cor all, of the polynomial coefficients, and the root analysis performed 
in terms of these narameters. This extension has proven exceptionally 
useful in two parameter analysis of control systems whose characteristic 
equations have coefficients that are each functions of the same two vari- 
ables or system parameters. 

The object of this thesis, then is to apply the results of Mitrovic 
and Siljak to the design of feedback compensation for the sampled-data 
control system. A secondary objective is to show the superiority of the 
two parameter method over the more classical techniques, and to demon- 
strate the adaptability of the entire design problem to solution by digital 


computer. 





eerie Komee Ons. 
Let us consider the characteristic equation of the sampled-data feed- 


back control system as, 


F(z) = 1 + GH(z) = 0 (la) 
able k 
or, Flz) = M a2 = 0 : (1b) 
=O 


where the a's are constant or parameter varying coefficients, and z is 
one of the complex z-plane locations of the closed loop roots of the system. 
Now, consider the generalized point in the z-plane, 
) 

nal Sh @ a: pozC 7 . 4 So 
where |z|= Ge 
- if 
cos@= j sing 
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If this general coordinate is substituted into equation (1b), the character- 


and arg z 


istic equation may then be written as: 


os K 
-» | 
F(z) = Ya,[o,§ +jOl- SFP & 
K=O z z 


Equation (3) is, in general, a non-zero quantity and assumes a zero value 
only when the value of dd, and iS assumed corresponds to that of one of 
the closed loop roots of the system. Thus, if we require that F(z) = 0, 
equation (3) describes the closed loop pole locations if the a's are 
constants, and every point on the z-plane root locus if one of the a's 


is parameter varying. Therefore, 
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2 ak 
a,_( WD, Ge oy (0, (1- S, )) 20 (4) 
K-38 


maps the characteristic equation from the z-plane to the Mitrovic plane if 


the 28 are constants, and the parameter plane if the a's are variable. 
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Only the parameter plane will be discussed hereafter. 


1,2 Conformal Mapping Techniques. 

The requirement that equation (4) be identically equal to zero 
implies that both the real and imaginary parts are zero. Siljak, [5], 
demonstrated that this magnitude requirement on equation (4) leads to 


the two simultaneous equations, 


te 
K 
— T, (Go aio (5a) 
and, 
i\ K = 
> & , CG) =o, (Sb) 


N= c. 


where T, ¢ oo) and Uf q. ) are the Chebyshev functions of the first and 


second kind respectively. Furthermore, since 


= Cc f a 
equations (5a) and (5b) can be rewritten as 


ey 


z. ~ ay WD, UL ee (7a) 
and oes 
_ ay an U.S) = 0 (7b) 


H Ry me aes 


respectively, with the UL ( Sere generated by the recursive Chebyshev 


formula 

Di eee UN >? ele ae, 6) 
with US? 20 

Ui ” - : 
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Siljak, [4] and [5], considered the case in which the a's are linear 


functions of two scalar parameters, 





a =Xb + Be +4, (9) 
The substitution of equation (9) into equations (7a) and (7b) forms the 


set of simultaneous equations, 


BIO, GS) Ee e.G) + (GS eno (10a) 


and 
HK BoC Way) + (Ss CoC 62,,6,) + Dg9(QzyO) = 9 » ~— (106) 
where 
K oe" K 
iia ea “bere se? “ 85s ¥ b, &, UG.) 
K=o0 K=O 
aN K rSu0 iK 
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Equations (10a) and (10b) can then be solved simultaneously for the two 


parameters, alpha and beta as; 


©,D, - C,D 5 B,D, - B,D 


2 o>, ~ BD, (11) 
en oan Be By) 








Equations (11) are the fundamental results of Siljak, [4] and [5], 
and provide the basis for the work which follows. Equation (11) can be 
interpreted as the calculated values of the two parameters alpha and beta 
required to put a closed loop root at a particularly chosen location on 
the z-plane. Thus, by fixing Os, and varying SS from -1 to +l, upper 
semi-circle contours of constant Ay can be mapped on the parameter plane. 
Thus, the argument of each upper semi-circle complex root conformally maps 
onto a single point on the parameter plane, the coordinates of which te 
the values of the two parameters producing that particular root configura- 
tion. Similiarly, in equation (11), S, can be fixed while G),is varied 
from zero to infinity producing loci of constant SC. » or angle, on the 


parameter plane. 





A third type of conformal mapping is yet required to completely specify 
all the roots of the characteristic equation on the parameter plane. The 
loci of constant GJ, and constant Su determine all the complex roots. 

The loci of points on the parameter plane corresponding to constant real 
roots on the z-plane can be obtained directly from the characteristic equa- 
tion, equation (lb). The substitution of equation (9) into equation (1b), 
and letting z= Gy , yields 
ww Q 
X y a+ BY ON + y 4 ee SS 
K20 K2O 


K=O 


which for a given J, is a straight line on the parameter plane, regard- 


less of the order of the characteristic equation. 





CHAPTER II 
SECOND ORDER SYSTEMS 

The second order sampled-data control system provides an excellent 
model for developing the feedback compensation techniques from applica- 
tions of the equations of Chapter I. For the initial design application, 
consider the pure second order system shown in figure (la). The selected 
design goal is to minimize the settling time by suitable choice of the feed- 
back parameters, a, and a,. If the states of the system of figure (la) are 


It 2 


defined as: 


ip 


X(t) c(t) 


il sad 


and Xy(t) = €(t) , (13) 
then the system can be represented by the signal flow-graph of figure (1b), 


and the vector-matrix equation 


X[(k+1)T) = @ xcr) + MReen), (14) 
where 
fe 2 
1 - (a, KT )/2 T - (a,kT )/2 
$ 7 - (a, KT) 1 - (a,KT) (15) 
and KT? /2 
ie = (16) 
KT 


The characteristic equation of the system is; 


Dee +- z(2a,KT + a kT” - 4) + (2 - 2a,KT Ss a kT’) = Q 


2.1 Constant Magnitude Loci. 


The duration of the transient response, or settling time for a 


linear system of the second order is, [6], 


a > § — 
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Figure 1 
Pure Second Order System 
Diagrams 





(a) Block Diagram 


é SS 





(bd) Signal Flow-Graph | 





t= 4/ = ee : (17) 


For a stable system then, the contours of lines of constant settling time 


of a second order system in the s-plane are 


c= 5.0, | 

or Ss = - 4/t ; (18) 
Accordingly, from one definition of the z-transformation, [2], 

z=e% ; (19) 


we can determine corresponding contours of constant settling time in the 
z-plane as, 
i Beaty. 


, (20) 
which, since T and e are constants for any particular design, are circles 
centered at the origin of the z-plane of radius ent coe 
Thus, the loci of constant Ws contours on the parameter plane can 

be interpreted as loci of constant settling times; and, when ()-2 and le 
are specified, the appropriate ae can be uniquely determined. Solving 
equation (20) for t > 

t. = - 4T/In Gs ‘ (21) 
Equation (21), however, must be interpreted in terms of the z-plane, dis- 
crete, characteristics. For example, for a step input, the minimum set- 
tling time obtainable is nT seconds, for an n-th order system, when all 
the closed loop roots are at the origin of the z-plane. Accordingly, 
when ts approaches zero in equation (21), t. should approach nT. Thus, 
for proper interpretation, equation (21) must be written as 

t= nt - 4T/In W, : (22) 
The particular conclusion of interest here is, in terms of the parameter 


plane, to minimize te one has only to minimize OJ, within the stable 
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region of the plane. 
To demonstrate the effectiveness of the constant W, loci in aid- 
ing the designer meet settling time criteria, consider again the pure 


second order system of figure (1), whose characteristic equation is given 


on page 6. It is first necessary to define the system parameters, so let 


2 
AK = akT (23) 


a ££ & ap KT 


By reference to equation (9) and the characteristic equation, then 


dy = 2 bo = 1 Com <2 
d, = -4 b, = ] Cc, = 2 
d., = 2 b, = 0 Co = 0 


which when substituted into equation (10) yields 
1s on Se ee OWS) 
oT ee ow oy gOS) 
Ae 
é 
ye WO, U, ( Ss COE. U, ¢ ~~ 


These factors can now be substituted into equation (11) and the constant 


\o) 
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ww) 
ti 


=, 
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Qe loci plotted with the aid of a desk calculator or a digital computer. 
Note that if a digital computer is used, as is most likely, a limiting 
process is required to remove the common factor of GOs. from the numera- 
tor and denominator of the parametric equations when attempting to locate 
the GQ), = O point on the parameter plane. Figure (2) shows the resulting 
lines of constant G)z on the parameter plane for this system. 

From figure (2), minimum settling time is predicted when i). = 0, 


or when 


1.0 , 


1) 


A. 
and Ve 


1.5 , 


1] —_ 
=_> —-™ 
= — 

li, 
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In order to determine the system response as predicted by equation (14), 


let 

T = K =9R0 : 
then a, = dy) 
and a5 = 1.5 is required. 


Table 1 lists the calculated magnitude values of the states of the system 
at consecutive sampling intervals using these values of the two parameters. 
The time response of the states is plotted on figure (3). A unit step in- 


put is assumed. 


-———_ = 
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Since in this example, n = 2, these results are in agreement with that 
predicted by equation (22), and the response is the minimum prototype, 
ripple free response. 

In order to demonstrate the most important design advantages in using 
the parameter plane for feedback compensation; assume that, for some other 
reason, one of the feedback coefficients is constant, say ay is fixed at 
2.25. Let K = T= 1.0, also be constant. With these constraints on the 
design problem, the parameter plane affords an excellent potential for 
setting the remaining degree of freedom to minimize the settling time, 
in so far as possible with feedback compensation alone. Referring to 
figure (2) for © = 2.25, Gz ,is minimum at 4), = 0.5 when = 1.875. 


From equation (22), 
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Figure 3 
Response of Pure Second Order syetan 
with Minimum Settling Time Compensation 
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is the predicted settling time. 
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Table 2 lists the state variable magni- 


tudes at the indicated sampling interval in response to a unit step input 


utilizing the above parameter values in equation (14). 


Figure (4) shows 


the time response of the state variables of the system as a function of 


the sampling interval, or time, 





TABLE 2 

i a ha yi 
0 0 Q 

} 0.5000 1.0000 
2 0.59000 -1.0000 
3 0.3750 0.7500 
A 0.5000 -0.5000 
5 0.4063 0.3125 
6 0.4688 -0,.1875 
7 0.4297 0.1094 
8 0.4531 -0,0625 
9 0.4395 0.0352 
10 0.4473 -0.0195 
1} 0.4429 0.0107 
LZ: 0.4453 -0,.0059 
is} 0.4440 0.0032 
14 0.4447 -0.0017 
iS 0.4443 0.0009 
16 0.4445 -0,0005 
er 0.4444 0.0003 
ies 0.4445 -0,0001 
19 0.4444 0.0001 
20 0.4444 -0,.0000 


The results indicated in Table 2 appear to discredit the validity of 


equation (22) in predicting settling time. Actually, equation (22) simply 
provides an engineering approximation to the settling time problem. Analy- 
sis of the results of Table 2 indicates the output position is within 1.95% 
of the final value after 8T seconds, although an actual zero velocity and 


final position of 0.4444 was not achieved until the 4lst sampling interval. 
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Figure h 
Response of Pure Second Order systen 
with Constraints 
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A similar test was performed with i 2.25 and = 1.50 for G), = 0.8. 
Eouation (22) predicts the settling time as t. = 20T, at which time the 
output had settled to within 3.3% of the final value. In that example, 
zero actual velocity did not occur until t = 99T. However, the settling 
time prediction of equation (22) remains within generally accepted good 


engineering design proximity of the desired value. It must be borne in 


mind however, that for a, different from 1.0, the conventional definition 


1 


of system error cannot be applied. Fortunately, in many cases, input/out- 
put scaling can be applied which permits the position voltage feedback to 


be set at other than unity. 


2,2 Constant Angle Loci. 


Equation (11), the parametric equations describing 


C_D - C.D 


19 2) , 
oe - Vix 
oT ee 2 


the Joci of the complex roots on the parameter plane are functions of both 


BD Bae. (11) 








B,C BAC, 
Ge and Ss . The constant magnitude curves were obtained by fixing 
Gas and allowing S to vary as the mapping parameter. Similiarly, 
Sx can be fixed while i), varies from zero to infinity to map lines of 
constant z-plane angle for the complex roots onto the parameter plane. 
While the loci of this constant angle on the parameter plane does not have 
any direct usefulness comparable to the loci of constant magnitude and the 
determination of settling time, they do lead to a very important secondary 
property. With loci of constant (J, and constant = drawn on the para- 
meter plane, every point on the z-plane now has a unique point defined on 


the parameter plane in the case of the second order systems. The fact that 
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for higher order systems, several z-plane points may map onto one para- 
meter plane point will be discussed later. 

For the second order case, this one-to-one correspondence between 
points on the z-plane and points on the parameter plane is useful in con- 
formal mapping of various loci from the z-plane to the parameter plane. 
Among the possibilities are the logarithmic constant damping spirals corres- 
ponding to lines of constant damping in the s-plane. Another use would be 
the mapping of lines of constant overshoot, [2] and [7], from the z-plane 
to the parameter plane. This z-plane mapping permits the designer to 
select feedback compensation which simultaneously satisfies overshoot and 


settling time criteria for the second order system. 


2.3 Constant Real Root Loci. 

In order to completely specify all possible z-plane root locations on 
the parameter plane, real root configurations must be allowed for as well 
as the complex root configurations. As indicated by equation (12), the 
loci of constant real roots are straight lines on the parameter plane. 

For the second order system, the closed loop poles are either both real or 
both complex. Therefore, the zone of real root description and the zone of 
complex root description cannot overlap, except at a critical point where 
the real root is about to emerge from the real axis. 

This zone separation characteristic greatly simplifies the sketching 
of the parameter plane for the second order system. Fortunately, only one 
set of calculations is necessary. The constant Gay curves must be cal- 
culated. However, since these loci of constant magnitude are straight 
lines, only a minimum amount of calculations is involved. In plotting 
the constant Ss, curves, one has only to connect the coordinates of 


equal S. on the ‘ °, curves. Then, to plot the lines of constant seen 
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since from equation (2) we can write, 
(eae IS | (24) 
one has only to draw straight lines tangent to the 1S | = 1.0 line at 


the OS, intersections. For example, a line drawn tangent to the S = -1.0 


line at the intersection of the CON = 0.5 line describes the real root 
locations Uz = -0.5. This sketching technique will not necessarily be 
extendable to the higher order systems unfortunately; however, it does 
increase the simplicity of the second order parameter plane. 

Figure (5) is the parameter plane for the pure second order system 
of figure 2 which now has some of the loci of constant angle (solid lines) 
and some of the loci of constant real roots (broken lines): superimposed 
upon it. Note that this parameter plane configuration also accurately 


describes the system stability in terms of the two parameters. 


2.4 Pole Effect. 

In order to study the pole effect of the second order systems on the 
parameter plane, the system shown in figure (6) was studied. This system 
can also be represented by the vector-matrix equation (14), where in this 


case 





| - SS (PTI +e") 2. a Kae(pT—1 +) 
7 @. Ka -PpT 


17 





oe oe 
sos 


t : = 5 Sao es t 
Siacete : eS eas a J a aS ts = — - fs = om (oaeied 5 i A ; ME EIBS, + te Fae an ane ee eS 
: Be 6 ER ee te a a Ye Sie ipa! Bruen GAN a: A 1b) bbeaba eps ele tes i Sa Tepes Lp aga Fe fal dig Aas ga by Sst esr AR re YES 
ie a | ee - —=_ ~~~ - - ~ —_ - -_ + + ee ee : > - 
o vabs ate $6 ~- aw bem wes 6 tro - = . 
cas eee pasereseneet 


' oe Pe Fe ne eg a ry : ’ 4 wate Hs 4 
. . omre : Dotlale ie lode bane oes r : 
ee a the ae u 
2 7 -~ 
: j--—1 pe 





= 





i i ' 2 Salts Ge ont ree tarlens = - . * ede tytn be 

s ; ~4-- 2 Se tee, oe oe oer ee 3 4 imei, te an ee td SS oe oore es ae eee - 

j ' po 5 . ca = 7 7 - ete} ee igccey tka a . ‘ > aaa ri + t= ~ J $5-2-.4- ‘ ~~ 
= <4 se . - - eo - 4 4 ed -4 = — 1 . -< “) oe -. - = 7-2 = } - - - a vom, 
s ~yw™ S as . 7 - - . - 4 - eee Ae, - 4 


‘ ~ a 
ore See te, om ne fie -eey ee eee : a). a a 


Se . ’ 3 . 
. ; ‘<a — i = = * =~ 
a = } - a . = _ =~ + - -~ - + .  —. ns ~ - ve 4 . -- 7 90 9- “ - 
. , - . ‘ + ’ if . ‘ Saw +7 ae Te 
- on- x -—F e in ae : po ae te ry ton inal de soctimeatingk Makeied : -- ene . , ‘ miches ot y = he - ‘ + 
2 ee lee ae > = ar : a5 - - - «. + , ee ae ee ee ee) = _- | ° efetwey te te = - 2 - Os 4--— « oo —— oo + ; - * u 
Stage JB) SSE eS een ee See ee ae, | || ere jem retcee phe 2a mee ——_ rot i i SO 
Suis ; - h y 4 = 
f e : ; ; : 2 | . ae ee 
. =-@e- *- - -_-e - 7. ~ 2-4 ’ . . es a A 
--q + ¢- . ee -i rw -e a Ts - - «= - a ths 


<= = : ~ Sf wee wand - =) ae : : = —— a carne 
rine wee de 3 s = = aE A SS -¢. + om t= gente gee 
: a5 8 ok Sif Sees ) Sanus eee 3 4 ae ae : * Saye is : is 
PE (ay ian 3 Coes aS w See é ree eae 











—~- 
—* 


4 -- ~ - + ‘«. - = a i 
- . ‘ Ly ‘ . 
+s -- -. ° — wee e . . =e 
. -~ $e . . : - -~ ~ , . ote eae | 
= =; <«4- - - fe -2 0 Ge “ - _ S 5 - in i a “ 2 a F ae > 
' 7 : < 
oe m ~-- -~ <= cee S eS ie - + - —_— : - - - ener 2) oq. el : m = 2 ere : 5 : = 
2a. | -7 ¥ ag , rr -~ aia 1 a Sed = ae . r . . r > lan 7 3 1 i +7 = = 4 ° : . : E . ‘ : ee) 
-- ; a | _—— riaeap wes ee oe o-—- : + --—- 4 = =e ies ee * + —4-= wy f=]. . + ‘ F 4 
. 4-4 2 ser quae = bane fo a - i. aha - ipo - “4 4 “$t + + iis : le-$- _ ‘ 7 — . 7 * . > 4 4 “— 
= - 4 ~ - = a ~ ~- . en - . — 5 te 3 se oe i oe - , t-< . 
° : . on = . : » 
: r : i ’ ‘ 
F ~ - : ---- — SSI = - eae 5 ; iia : * . 
: eg {is : ‘ : - 
= 7 . ore -+ «— o-+- = . @ tee 4 
T : 






a ee 
ee S3pse reese a 
> OLS ALAS Baa paisaty Ysviatea es 


Es 4 4ory eee eee 
oy, ~—t = - 


—- 





el ee eee ee 


~~. yt me ps i ttt te ee tei -l(CQQlCU eee ee Oe le De 
aapemerge Os Me Bc ae 


++ + 


Feces gece: bescs tasassassvesest 
SES Eset 


aie aaa ie rstrttr4 suaneee - 


DAMCSASOCwH9AseAe Leg 
1 kAe POSES Peebe ft Fae PEs 
< 


i - a senthins at. - 5 eer pos cesEE APACE SSeHsT — 
—s rt + + + ) 4-4 bres os - = : _| Leen S208 as . 
: eee {3 , ¥ ae : : a a = ‘ < apane 

: 5 S mn 
. t . - = one ~ 
- a os « a ° - A, - e oe = = 
Y = 


eo ae 


- ~~ 


ew 


te te 4 _ : _- - - 4: -- t<--y - Soe Peale 
2 =i : SS eee tee a nd oad - boo = = Aote . F inbeg i, ee ee Stas == n = = A 2 : . oes 
wae tot eet ap rary tm + peers Fae a a 5 mee : = saa a eae ane ie a ! . : 53 es ate = me ie . ~ - = - Sa -- ee a 
TUS A} (Sa CG Pa Me foe Raped a es CS tian Boat: i Ly eeieare 3 oy : :: pai " ee : " i ; : TT epsn 
fies cs ) : me : 2 op LY 
volar ate ee lee ¢ pe an i i is * 5 : a F ree Ve : ri. betas = : . z -%, 
~4 op ee - - — ie - - -“——- + E 
= nn eo ———~ ~ = 5 . < Sees ae ee -$ owen. = 


: ote a % 
ih a 
eee a -_=< T°" 4 - — 9 > 
ees 


ae TEP 


* 
| neat trial oie oe 
u rf ; ' iF 





el Pee eee 
die eS ae [a 


. ee 
= 5 oe wl Sw yer ~~~ - Sos = 3 Z =a . ’ ne a a : rom = = ° a = " ; 5 

as oes a earner es | eee eee ke ee a = a - a t 

-- = +t -% ent OT mee ee eth et ere ee - oa : 


+--+ Ls ee bebmany pate tom oed men 


a4 444 (oo — i 


se fee Bite bet Tod 7 7 








- 
+ $= 2 varee ebb owe 5 a i 
+ & ~+ tee 
at - —pryze BET Sd a 
re — +t 
t t+ Ftacest Fad 
a 
-4- oqo bp, bdo ie. wee ~ -; : om yep - “4 - an | eS Sa, 
+ eo =—-f-eee = ww gue bee wre a ee is — + noe ~ 2 re £ . --e oe ee ~~ =: . - 4 t ie ~ . Fen fore lh foros -— = 
_— es bs) 6 awe - =o. - re ee arses a ray me = { wale — A So wee Oe : - es (eos ts - 
ore - @8- 9 oee-- +> sdbdichs Tntntn catia aie -- : - er = ae ee ee 4 +- pore - be yereny - 
+ =a - © ne - ee ee = =e -~ A nn ° cm ee - -—- - “, 2 ore 
5 = Siete ean a. . —s 2 _ 
5 £ + @be = 4 +~ r 
- ’ : = - = : +° -fen-e eters fin = ’ : - d o be ton fee 4 Z Fs 
‘ = ; ; oo, 
2 : 4 535 = fs e- - 5 5 
ee te ee ee em ee ee ee —_—— . 
t oof , 





- — 2 a - ee oo - | 
- 4 J c ’ 5 : x, - 
' RS +] ; - * : ; = ES oa -- mi — “ 4 = -* eo det ha I 
3 id d i a . Hl ine s i : " uy a : = as ry = ehotoeas : ip eee birt) (Res «. 
wm ee . po cee Se oh . ne Se See ee a Oe , ee eee Sa Ase ~)_ ee cos o HRs SS eS 
} -. . - wm . y n ee ah ot O 2 2” - ‘ ~$ « + = aie —-e- e . - 24 . 4- es - oy t edn a 
= mee 4). ha ay eae f . aa A Z es a a -° a-eo@ee — - esis Soe - emt wee Fo ee py a, 
ee od - 8 . whe aoe =f ise tee ec ead we + @enecasceb-- - Ck one . im 1 4-- - -4de.. “0-4 - ~ r - «e4-=+P mpm dye 5 
} - ; e ie BS come - a G a : 3 SH he, 2 t- ‘| Pom’ z : ae ee ‘4 = Pie s : ee 5 at al ea ee Ee Pies ices 
Ce ener cor ear meen ee ere ne re ee qh ee re ee ee Se SS ee en ay ER Oh ee ee ee i er a SO Ot we - eee py ee he ee eto 











Figure 6 
General Second Order 
System Diagrams 





(a) Block Diagram 
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(bd) Signal Flow-Graph 
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The characteristic equation is 


2 . 
z t+ A,2 + Ay = 0 (275) 


where 
A, = Ka, (T/p-1/pte™/p) + Ka,(1/p-e"P?/p)-1-e P 


and 
Ay = Ka, (1/p>-Te P"/p-e P*/p*) + Ka,(eP*/p-1/p)teP*. 
One important difference between the pure second order system and the 
system with damping becomes evident by inspecting the characteristic equa- 
tion above. In the pure second order case, by judicious definition of the 
system parameters, equation (23), it was possible to produce a parameter 
plane as a function (indirectly) of all physical parameters, K, T, ay» and 


a By moving one of the poles into the left-half of the s-plane, it is 


9° 
no longer possible to obtain parameters involving T. However, the system 
gain, K, remains a common factor in the characteristic equation coeffici- 
ents. Therefore, for the second order system with damping, we can define 


the parameters as 


ee Ka, (26) 
yo - Kay 


For further study of the pole effects, let p = T= 1.0, then from equa- 


i> 


tion (9), 
dy = (),368 by = 0,264 co = -0,632 
d, = -], 368 by = 0.368 c, = 0.632 
d. = 1.0 b, = 0 Co = ( 





and from equations (10), 


B, = -0.264U_,( S_) By = 0.368 6), 0, ¢ ) 
Cc, = 0.6320 (SS ) Cy = 0.632 GU, ( S. ) 


] 
tt 


4 
ie co ccu SS, “HO, Y, « Ss? 
-1,368 QJ, ( S) + CN" U, ( 5). 


o 
tt 


Figure (7) is the parameter plane for this system, with loci plotted 
for comparison with figure (5). Note that the shifting of one pole from 
the origin of the s-plane into the left-half plane has produced a counter- 
clockwise rotation and a compression of the upper semi-circle contours. 
This suggests that the further into the left half plane of the real pole 
(closer to the origin of the z-plane) the more sensitive the system re- 
sponse is to changes in ay Additionally, from figure (7), the extremities 
of the complex root region on the parameter plane have been extended, to 
the point where negative values of beta are now within the stable region. 
From figure (7), minimum prototype response to a step input should be ex- 


pected for 


HA = 1.58 
iS 


Consequently, there are two ways to achieve minimum prototype compensa- 


1.24 


tion. First, if K is fixed to some value other than 1.58, for example 
1.0, then non-unity position feedback must be used. With K = 1.0, the 
magnitudes of the state variables at consecutive sampling intervals in 
response to a unit step input, calculated using equation (14), are shown 
by table 3. Note that this response is identical to that shown on figure 
(3), except for a constant scale factor, which is naturally equal to 


1/a,. 
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TABLE 3 


nt X(t) xX, (t) 








0 

0.368 
Om 32 
0.632 
0.632 
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For the particular case where the conventional error definition is 
required, to obtain the minimum prototype response, K must be restricted 


to the value of 1.58. Hence, the corresponding value for a, is 1.24/K or 


2 
0.785. The calculated response using these parameter settings is shown 
by table 4, assuming a unit step input again. 


TABLE 4 


nT X, (t) X, (t) 








0 

0.582 
1.000 
1.000 
1.000 


EWN KH OS 
Oo Oo 62 = © 


2.5 Zero Effect. 
To study the effect of a real zero on the second order parameter 
plane, the system shown in figure (8) was analyzed. For this system, 


the state transition matrices are: 


1 - Ka, (T + 27° /2) T - Ka, (T + 20° /2) 


= Ka, G1 + ZT) 1 - Ka, (1 + ZT) 
and 
2 
K(T + ZT /2) 
K(1 + ZT) 
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Figure 8 
Second Order System 
with Zero, 
Diagrams 
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Again, the characteristic equation is 


2 
Z + Aj2 + Ao = 0 : 
where 
Ay = Ka, (21° - 2T) - Ka, (22T +2) +1 
2 
and A, = Ka, (ZT + 2T) + Ka, (22T +2) - 2 - 


Again, a choice of parameters for the parameter plane can be made as was 
performed in section 2.4. For further study, Z and T must be fixed in 
magnitude. For this example then, let Z= T = 1.0. Following the now 
familiar calculations, the parameter plane of figure (9) results. Similiar 
to the results shown for the pole effect, the parameter plane loci have 
been rotated counter-clockwise from the pure second order configuration of 
figure (5). Apparently this counter-clockwise rotation can be associated 
with the relative stability of the system, as both the system with the zero 
and the system with the left-half s-plane pole are inherently more stable 
than the pure second order servo. However, no measure for quantitatively 
describing this relative stability has been discovered. 

From figure (9), minimum prototype, ripple free response to a step in- 


put is predicted for 


if 


A 
fp 


Again there are two choices for the physical parameters in achieving 


0.5 


OF 125 


this compensation. If K is variable, unity position feedback may be used. 
If K is fixed, then non-unity position feedback is required. Responses 
similiar to that achieved in the last section result for each of these 


options. 
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CHAPTER IIT 


HIGHER ORDER SYSTEMS 


3.1 Parameter Limitations. 

It can be readily demonstrated that in order to achieve complete free- 
dom in locating the closed loop roots of an n-th order characteristic 
equation, the control system error signal must contain feedback intelli- 
gence on each of the n dynamically independent "states" of the system. 
Consequently, the system analog or mathematical model consists of n para- 
meters and the parameter space becomes an n-space. Unfortunately, the 
existing Mitrovic, or algebraic, methods have their basis in the solution 
of two simultaneous equations, equations (5a) and (5b). This restricts 
the mathematical treatment of any n-th order system to only two degrees of 
freedom. While this reduction in the degrees of freedom restricts the 
ability to position the closed loop poles in the z-plane, a large class 
of control systems can be reduced to the two parameter problem. 

Typically, the gain setting, K, in figure (10) is dictated by the 
requirements for steady state error. A second degree of freedom is normal- 
ly eliminated by the adoption of the conventional error definition which 
demands unity feedback. Finally, in most physical systems, noise free 
derivative signals higher than the second derivative are unobtainable. 
Thus, the only degrees of freedom for feedback compensation typically avail- 
able are the magnitudes of the first and second derivative feedback. 

With this restriction on the positioning of the closed loop roots in 
mind then, the feedback compensation of systems with order greater than two 
can be designed on the parameter plane using techniques already discussed. 


Two additional features of the parameter plane must now be considered. 


ZT 





Figure 10 
Typical Block Diagram for 
n-th Order Sampled-Data System 
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Unfortunately, at the present time, these two features tend to reduce 
the effectiveness of the Mitrovic Method. However, it is believed that 
this reduced effectiveness is still an improvement over the existing one 


parameter design techniques. 


3.2 Maximum Modulus Uncertainity. 

The parametric equations describing the parameter plane loci were de- 
rived from the system characteristic equation. Any specified location (M- 
point) on the resulting parameter plane can be interpreted as the values 
of the two parameters which will put a closed loop root at the specified 
position of the z-plane. Therefore, for a typical fourth order system 
for example, the closed loop root configurations shown by figures(lla) and 
(11b} each correspond to a unique point on the parameter plane, each lying 
on the Od. = 1,0 loci. In order to distinguish between the different con- 
figurations on the parameter plane, it is also necessary to plot the = O2D 
and i), = 1.5 loci. The intersection of two of the loci of constant ive 
then adequately describes the location of each pair of complex roots. If 
figure (12) represents the parameter plane for the fourth order system of 
figure (11), then the root location depicted by figure (lla) corresponds 
to point A, while the configuration of figure (11b) corresponds to point 
B on the parameter plane. 

Thus, for an n-th order system, the modulus of each real root or each 
pair of complex roots is completely determined on the parameter plane. 

The accuracy of that determination, however, is directly related to the 
nunber of loci of constant OJ, drawn. In the second order case, it was 
fairly easy to extrapolate between adjacent loci. As the order of the 
characteristic equation increases, the difficulty of making this extra- 


polation increases. It should be noted that this problem of interpreting 
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Figure ll 
Tenia z-Plane Closed Loop Root 
Configuration for 
Fourth Order System 


Figure (lla) 





Figure (11b) : 





Figure 12 
Typical Parameter Plane 
for Fourth Order System 
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intersections of constant loci also exists in the second order case when 
the plant gain is insufficient to produce complex roots. Uncertainity 
arises as the order of the system increases. For example, given a fifth 
order system, and a corresponding parameter plane M-point defined by the 
intersection of two constant i), loci and one real root loci. Obviously, 
this M-point must determine two pair of complex roots and one real root. 


However, with existing techniques, it is impossible to associate any 


specific c), with a specific complex root. 


eo. Discontinuities on the Parameter Plane. 
A second feature which may exist in higher order systems involves dis- 
continuities on the constant ‘.', and . loci. Equation (11) describes 


the parametric equations for alpha and beta as a ratio of two polynomials. 


Consequently, whenever 


BLS, 


= = 0 6 
BoC, d (26) 

a discontinuity exists on the loci of constant 6), or S . At the present 

time, the physical significance of this discontinuity is not understood, 


in fact one may not even exist. Its presence, however, does complicate 


the sketching and interpretation of the parameter plane. 


3.4 Dominance and Settling Time. 


In the case of the linear, continuous system, the response of a 
higher order system can often be discussed in terms of a dominant pair of 
roots in the s-plane. This is a pair of complex conjugate roots whose 
residue and/or transient settling time leads to by far the greatest 
portion of the observed physical response. In the s-plane, these roots 
have an argument and a real part considerably less in magnitude than any 


other pair of complex, or real, roots. If such a pair of roots exist, the 


a2 





response can be quite accurately predicted in terms of a second order 
approximation and the settling time predicted by equation (17). 

Similiarly, for the sampled-data system, the dominance of the system 
response can be discussed in terms of the pair of complex roots, or real 
root, whose modulus or argument is conside: ably greater than that of any 
other root, real or complex. Jf such a root exists, the settling time 
can be closely approximated by a slight modification to equation (22) such 
that 

Cos nT - 4T/1n ( ew (27) 
When the dominance cannot be clearly determined by the difference in the 
various arguments, equation (27) becomes less and less accurate in pre- 
dicting the system settline time. 

Alternatively, the minimization of the settling time can be achieved 
qualitatively by minimizing the modulus of each root or pair of complex 
roots. That is to say, compensation must be chosen which will place all 
the z-plane roots within a circle of minimum radius. If this circle of 
minimum radius happens to be the origin of the z-plane, then minimum pro- 
totype response is guaranteed. If, however, all the z-plane roots cannot 
be placed at the origin due to the limitations on the degrees of freedom 
available, minimum settling time will be somewhat greater than minimum 
prototype, and will not be ripple free. On the parameter plane, this 
minimization of the settling time involves the examination of the inter- 
sections of the constant bo, loci for each permissible M-point, and 
selecting that one which does place all the roots within the circle of 


minimum possible radius. Observations of the resulting arguments will 


enable the designer to determine if equation (27) is then valid or not. 


55 





235.) lbied Order Example. 


Consider the system shown on figure (13). This system can be re- 


presented by the signal flow graph of figure (14), and the vector-matrix 


equation 
x fa+it) = f xr + A eck) 
where 
Py 3/2 - 2077 + 12e727) Casa - 7 + 1/2072) 
$= 0 (em - oN (a - a ou 
0 (De 3 + ye ae (<2@ + deb) 


T/2 - 3/4 + at - Wace 


L722 a i/a(er} 


lc 


and, e(kT) R(kT) -X, (kT) -a,X, (kT) - aX, (kT) 


The characteristic equation, evaluated at T = 1 second is 


3 2 

Z + Aoz +A 2+ A, = 0 : 
with A, = 0.1998a, + 0.2325a, - 1.4191 

A, = -0.1263a, - 0.4650a,, + 0.7238 


and A -0.0734a, + 0.2325a, - 0.0310 ; 


0 i 

Observe that it is no longer possible, as was done for the second order 
systems, to include the plant gain, K, in the definition of the system 
parameters. This combination of the gain, K, with the feedback coeffici- 
ents is only possible so long as each coefficient in the numerator of the 
product G(s)H(s) is some product containing a variable feedback gain. In 


the second order case, 


G(s)H(s) = K(a, + ays) 
D(s) 





The characteristic equation can be formed as 
1 + z-transform of FG(s)H(s) | 
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Thus, the only place that the feedback coefficients occur in the character- 
istic equation, they are each multiplied by the gain, K. However, for 


this third order example, the produce G(s)H(s) has the form, 


G(s)H(s) = KC 1 + as + a,s-) | 
D (s) 
In the characteristic equation then, K is a multiplier of terms other than 
those involving the feedback coefficients, ay and Ay - For the example, we 


are restricted to defining 


i> 


a 


A 1 
2 ay) 


Following the usual manipulations, the third order parameter plane, figure 


lif 


(15), was achieved. No lines of constant S, Were included on the para- 
meter plane, only to avoid unnecessary cluttering. Observe that the separa- 
tion of complex root regions and real root regions no longer exists; which 
should be expected, as the root locus always has one real root. A/so note 
that the region of stability is that region of the parameter plane which 
is enclosed by the contours C= 1.0, Op = -1.0, IES = 1.0. From figure 
(15), minimum settling time using feedback compensation is predicted when 
a, = 1.61 
a5 = 0.54 
The circle of minimum radius enclosing the roots is G), = O23; 
Table 5 lists the calculated magnitudes of the state variables at consecu- 
tive sampling intervals in response to a unit step input using these values 
for the feedback coefficients. From Equation (27) then, minimum settling 


time is calculated as 


3T - 4T/in 0.3 


ct 
tt 


= 6.3T 


37 











te ion ahis Eenard + Ee - irr 
: erssoeeswesss Rape, ye eae 














yoy ye 


= i 
ee Se ee 
a ; i) lee tel 2 fe 
~ - olen o> Sunes - a4 4e-)- 





“e oe eon 
~~ er = - 








ee : eee ere = 


we Pe a 








SE bewossioos fare 


Rene re ht he eT 





= = - 
Set met oh ed + . 
Boe cae. 7 agaee 
: w§~s4 . A “ : 
aot —~ « ad 1 . + 2 ; 
— ‘ 

. - oie 

. & 

= ~“ 4 


nes 


Sead secetosst tease 
shane SeeRy &@aos 








lasses . 
-i- igs 
, i-- . ei 
a . > 5 a= - 2 
~~ — we top er 0 fm —— « we hee 
pw wo er ee eee . - tes 
aioe 
- + - — - - : 
‘ 
2 - © abes ba =e 
° « 


eS) Ea 
oe Regis cans ° 
t+ 


sean 














i Son9 pees ion owes ++-- 


See 
Eire § Pep +43 SES 


e- 
u R at Be Sess Soeras 
a ft Om Real | 4 a : a. > 5 ae. ae et t 
eedlenses : r. 2 ag "Tag  spebe Boo sasus genus : 
ae nee Ty a ere 5 $ ~~ Ji ose aa dt j-te ry ro ey ere | ty fates hee ee ft ae 
> @ ~~ .- bra . , ° - . - ' © geet bee a = ta SB day wtp Fe poe ty att te tg a ated ee ae ad be - 
: = eee SS Saree : Si oe 
- - . “ mat 5 
. z . E 












comes ba - _ ap ane * r 
hie | Pea AW (ot 2o ee a ee pipe Bl Rae A ee dies LA «Lap me wht ejtaon aera 
“> Fey aoe pnb ada: 4 “ew oben, o- mt jose sens 1 hg fe te th nape fe ae F- sana Se Be re = 
5 j:=- ote ~~ —t apa * - + reo i bam smth fin lamin 8 = +e war Ong jenn me t=? oop) bir — — a=: 
= mare Be oe a pS ae Ee PEE. * rtos aah hese 4 Jp bata mn 7 —-+ 


Sot i “aT, rere ee | = & Sie ran s Fs 2 + [ eter tie ~~ te: sree eee 4-9 ie = 4 5-6 = — ae a ery a eo 
“1 BS) eka “4 atlas | t 3 ae - 4° -4 bs se 44 tes = tr re os, ee | Bree 4 - i - Oe. Ite ci oo ae. 
—~- } a - boot eel E 4++ +44 + ie a 






























; x frie t 1 
« og J. 4d -+-9— 6 ie — “+ 3 + « eet dew og bf apad ng > 
Siege Sha 1 ee tt 7: i= poles eles es 4 as fe 
7 fiisgce: a pe ene bee pr antaaena ‘oy—then : ae ne os —s egies - 
— Pe | - - | peg oes? aes OS “p> + 404 oe 
oe ee es io x 4-4 (s- 
~ + > 


{ vo ae | ae + +: 
°° - +. + ia -p= { . 








pea bes pas bu ee 





‘ i 7 fis a. wt 7° 2 
+ ' =< . 
- 4-4 * i tee (esi 
. . - }- -+> 4. + +4 4 ;2 ae 
= oe . * a tt ' Bs 
+ 4- t . 2 {- 
+ 





- + > 4 = - . 
rf _- res - 4a oe 
aes eee <b an ee ee ee ody be aa 
i; - jam tat pes eee sera > 
ec: i 7 en mee _ 

















= =e bday 
os ellie a crn Bee 
| ee 8: fe . 
a on 2 eS ee eee a ee ee 
=} _——— 4 -+—-- 
aaa | —r-¢ i -=t 942 eee Tet ee 
oe Wyo =: a ged ge 
a | 3 _—- ~+h ew —- 
SR slams 
ee iii ot 
©] - | : eee 
~ VW 
ee ge eee ee 





OT 3G By ait 1a i ora 5 are a) ie m7 A ot | OP Aas ary ; VWyapeditis iv ittt S225) ae © 
eee ime za re aman i Mina detec wari pe ee Ae peed ro Sst) eae 


6 
arte ee * - 
T° t- ermigraanaesy ts yt on? + oS ue = oo xe -- _ y bee pew le + - 
+. = Soar peewee eee + gate e+ —— es ) SlusEs St = awe <j Be el a ee oe —- re aaa <5 ene pels 
esi 6 - Oe bee i wee “ et ane T wy wed wet ot 4 ome a= ee. -_4 = qed” obi id - aajee wate o—t ¢- 
=] - 44 = 4-6 Pr é . 
Bs . = . +4 = 7 ~- 4 wot oon ~4e4-d ~< 
-_- see “« 
- ~— ae) -~ — 
“; a g t = == ,2e~ 
2 ee ee ° 
A a eT ¥ 3 : 
SoS ee ae: 
hey gwen ee ey a Se eee —. - ae Y 
4 ‘ane ae eee b+ eae = 4: se 
ean jee f Jeewalf> - 









At time 7T, table 5 shows that the output has settled to within 5% of the 
final value. Figure 16 shows the state variable time response for this 
minimum settling time compensation. These results verify the effective- 
ness of equation (27) in predicting the servo settling time to a good 
engineering approximation, even with the dominance of the roots at maxi- 


mum modulus, as discussed in section 3.4, somewhat uncertain. 











TABLE 5 

aE X, (t) X(t) X(t) 

0 0 0 0 

1 0.08405 0.19979 0.23254 
2 0.33610 0.26768 -0.00653 
3 0.57736 0.20644 -0.06887 
4 0.74605 0.13340 -0.05968 
5 0.85111 0.08048 -0.03964 
6 0.91357 0.04724 -0.02411 
7 0.95002 0.02743 -0.01419 
8 0.97114 0.01586 -0.00824 
9 0.98334 0.00916 -0.00477 
10 0.99039 0.00529 -0.00275 
11 0.99445 0.00305 -0.00159 
12 0.99680 0.00176 -0.00092 
13 0.99815 0.00102 -0.00053 
14 0.99893 0.00059 -0.00031 
15 0.99938 0.00034 -0.00018 
16 0.99965 0.00020 -0.00010 
17 0.99980 0.00011 -0.00006 
18 0.99988 0.00007 -0.00003 
19 0.99993 0.00004 -0.00002 
20 0.99999 0.00001 -0.00001 


550 Algebraic Design Methods. 


One simple and quite obvious method for designing feedback compensa- 
tion for minimum prototype response should not be overlooked. In determin- 
ing the system characteristic equation as a function of the feedback vari- 
ables, which must be done to apply the Mitrovic Method, an alternate method 
arises for determining minimum prototype compensation. To obtain minimum 


prototype response for a given system, the function of any compensation 
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Figure 16 
State Variable Time Response for Parameter 
Plane Compensated Third Order System 
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must be to drive all the closed loop roots to the origin of the z-plane. 
The characteristic equation of an n-th order system, such as shown on 


Pi cure.G) 0). 1s 


Zz +A {2 Tr eo @?@e@6@6h6UhhrhUC< MMC NHhC<“ CO + Ay = 0 (28) 
where 

A. = Boi 4 + Bi54) aD A ae er Bien + BS 
and a. = feedback coefficient of the i-th derivative 


j 
feedback. Therefore, to algebraically design the minimum prototype compen- 


sation, thereby positioning each closed loop root at the origin of the z- 


plane, one has only to simultaneously solve the n non-homogeneous equa- 


tions 
Al = 0 
A129 = 0 
Ay = 0 


The solution to these simultaneous equations determines the values of the 
n feedback parameters required to produce minimum prototype response. 
Additionally, as was pointed out in section 3.5, since the number of feed- 
back variables now equals the degree of the characteristic equation, we 
can rewrite the A. in equation (28) as 

A. = Bsa + B, Ka, gos 5) et By ka, + BS F 
thereby providing the designer with yet another potential variable. 

To demonstrate the use of this algebraic design technique, consider 
again the third order example of figure (13). Let the open loop gain be 
variable, K, instead of 1.0 as shown. Furthermore, let the coefficient of 
the position feedback path be a,. Then, for minimum prototype response, 


0 


the characteristic equation becomes 
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e + A 22 +A.z2+ A. = O 


2 1 0 
where A, = 0.2325Ka, sip 0.1998Ka, 1% 0.0840Ka, - 1.503 = 0 
A, = -0.4650Ka, - 0.1263Ka, + 0.1709Ka. + 0.5529 = 0 


1 Z 1 0 


0 = 0.2325Ka, = 0.0734Ka, + 0.0188Ka, - 0.0498 = 0 


> 
| 


The simultaneous solution of these equations then, is 


Kay = 3,652 
Ka, = 4,448 
Ka, = 1.323 


In order to permit the adoption of the conventional error definition, 


then, KK =7 32652 
ay = 1.0 
ay = 1.218 
ay = 0.363 


Using these values of the system parameters, the response of the system 
of figure (10) to a unit step input is as shown by table 6. The time 


response of the state variables is shown on figure (17). 











TABLE 6 
a XK, (t) X(t) X(t) 
0 0.000 0.000 0.000 
1 0.307 0.730 0.849 
2 0.929 0.286 -0.849 
3 1.000 0.000 0.000 
4 1.000 0.000 0.000 
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Figure 17 
State Variable Time Response for Algebraic 
Compensated Third Order System 
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CHAPTER IV 
DIGITAL COMPUTATION IN MULTI-ORDERED PARAMETER SPACE 

While the feedback compensation of many n-th order sampled-data con- 
trol systems can be reduced to the two parameter problem presented in 
Chapter III, circumstances can, and do, arise where derivative signals 
higher than the second derivative are available. Since each additional 
degree of freedom available to the designer enlarges the region in the z- 
plane in which the closed loop roots may be positioned, common-sense alone 
dictates their usage where available. As pointed out in preceding chapters, 
Mitrovic's Method is primarily two parameter in nature. Additional design 
techniques are then desirable whereby a parameter space higher than a 2- 
space can be achieved. 

The digital computer provides one method of increasing the order of 
the parameter space. Any design process or technique is simply a method 
for locating the roots of a polynomial, the system characteristic equation, 
as a function of one or more parameters. These methods were developed pri- 
marily to surmount the difficulty of factoring a given polynomial. With 
the digital computer however, this polynomial factoring can be repeatedly 
performed, always with greater accuracy, and sometimes in less total time 
than analysis performed by the graphical and semi-graphical techniques. In 
this chapter, one approach to multi-parameter study using the digital com- 
puter will be presented. The emphasis on the application of this method 
will be directed toward feedback compensation; however, the extension to 
cascade compensation or other parameter variations does exist. 

Recent emphasis on the state space method of analysis, [2], as applied 


to sampled-data systems affords an excellent, if not necessary, point of 
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departure for digital computer analysis. The general form of the state 
space vector-matrix equation describing a sampled-data feedback control 
system was given by equation (14), 
X (e+ DT) = B xcer) + A R(KT) 

It can easily be demonstrated, see appendix 1, that the characteristic 
equation of the state transition matrix, g, is indeed the characteristic 
equation of the closed-loop transfer function. Since numerous numerical 
methods presently exist for determining the eigenvalues of a given matrix, 
the determination of the closed loop pole locations by evaluating the 
eigenvalues of the $ matrix is readily performed on the digital computer. 
To determine these root locations as a function of one or more parameters, 
one has only to repeatedly increment the parameters and determine the eigen- 
values of the resulting f matrix. In this fashion, an n-coordinate 
space can readily be created by computer programming techniques. One 
method is to use the FORTRAN language, and program DO LOOP's inside other 
DO LOOP's until sufficient coordinates are defined to permit the designer 
to sweep through a desired region in the n-space, determining the root 
locations corresponding to each point in the n dimensional parameter space. 

In order to demonstrate this method of analysis more fully, the 
following specific project was undertaken. Given the perfectly general 
second order sampled-data system shown in figure (18), the values of the 
varameters Ka, and Ka, required for minimum prototype, minimum settling time 
response to a step input are to be determined as functions of the sampling 
interval, T. The resulting computer program should be such that the de- 
sired design data can be readily generated with only a minimum amount of 
additional input data supplied to the data, for any second order system 
with real poles. 
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The system of figure (18) can be represented by the open-loop vector 
matrix equation: 
X[(k + 1)T) = FX(kT) + A e(kT). 
The error signal at the k-th sampling interval is represented as; 
e(kT) = R(kT) - a,X, (kT) - aX, (kT) 
In the chosen example, the F matrix is perfectly general, irregardless of 


the location of the plant poles and zeros, and 


p e Py! -p abo. ani - ene 
2 1 
F = % 

(pote paye eel cone 

9, Ce 1 2. ~ 1)/% P,e 1” + Pre 2 
where 
Say 
% oon: (P5-P,) 


A 
a) = 1, D- or 

Modifications to the F matrix are needed, however, if a pure second order 

system is used. The matrix on the other hand, is dependent upon the 

number of pure integration: in the forward path of the plant transfer 


Function. For example, for pn = Os 


ae 
= ; 2 
b, + b,T l. bot /2 
B =k 


»§@ ++ b, 4. bot 


When only one pure integration exists, that is P, 0, p, = 0, 


2. P -p,T 
hs (1/p,)° [(b, + bpT)p, = by + 90°F 1°] 
AA -p T 
A : P 
b, V(r) + B /p, > Qee layed 
where, 
0, = (b b.)p, +b 
49 POP) oa i 
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Figure 18 
General Second Order 
System With Zeros . 
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The program which was developed, see appendix 2, to satisfy the 
design goal requires four inputs in the form of data. In the order in 
which they are required, they are: 

a. the pole locations of the plant in the s-plane. The nega- 
tive value of this s-plane location must be entered on the data card, and 
the poles should be listed in decreasing magnitude order. 

b. the plant s-plane numerator coefficients listed in ascend- 
ing order. 

c. the minimum and maximum values of T, the sampling interval, 
and the incremental value of the sampling interval to be used in the in- 
vestigation. T minimum must be greater than zero. 

d. miscellaneous title information and graph scales required 
for graph output in accordance with an existing graph subroutine (DRAW). 

The program outputs are: 

a. a tabulated listing of the minimum settling time values of 
Ka, and Ka, determined for each value of T used in the investigation. 


1 


b. a graph plot of the minimum settling time value of Ka, vs. 


c. a graph plot of the minimum settling time value of Ka, vs. 
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Let us now consider the geometric significance of this program. It 
has been demonstrated that for a given value of the sampling interval T, 
the origin of the z-plane maps into a single point on the parameter plane 
for a second order system. However, if a three dimensional right-handed 
coordinate set is created, whose axes are Ka,» Kay, and T, the minimum 
settling time point-corresponding to CD, = 0 will map into a curve in the 
3-space. The plots of Ka, and Ka, required for minimum settling time re- 
sponse vs. T can be thought of as the projection of this curve in 3-space 
onto the Ka._-T and a) planes. Similiarly, any n-th order system, with 


HI 


m degrees of freedom, where m < n, can be investigated in a computer 
oriented m-space and displayed to the designer by a set of m projections 
on the a. VS. T, (i = 1, 2,...,m), planes. If m =n, a set of projections 
on the Ka, -T planes, (4°-= 1, 2:2.-,0), are available: 

The following examples will serve to illustrate the class of feed- 
back compensation problems, the solutions of which can be greatly simpli- 
fied by using this particular programming technique. 

EXAMPLE #1 
A sampled-data control system whose plant transfer function is 
G(s) = K(s + 1)/s(s + 2) : 
is to be used with a dual capacity for sampling rates of 10cps and leps. 
Compensation is to be designed to provide minimum prototype response to 
a step input for each of the possible sampling frequencies. 

Figures (19) and (20) are the minimum settling time parameter planes 
for this plant transfer function. From these parameter planes, the re- 
quired parameter values for each sampling frequency can be obtained as: 

Ka Ka 


10.02 0.949 
ims 0.538 


2 
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Finally, for unity position feedback, the physical plant parameters 


must be: 
i Ka Ka, 
0.1 sec 10.0 0.095 
1.0 sec | Pe ee 0.467 


Figure (21) represents the block diagram of the compensated system. 
When the three switches in the physical plant are "ganged", then the control- 
led switch of the sampling frequency will introduce the correct, associated 


feedback compensation for minimum settling time. 
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EXAMPLE #2 
A system whose plant transfer function is 
G(s) = K(s + Z)/s(s + 2) 

is to be controlled by a sampled error signal. The plant zero location 
is a function of a cascaded compensator, and can be placed anywhere on 
the real axis of the s-plane. A sampling interval of 0.5 seconds is 
intended. The designer is to set the gain, K, the zero location, Z, and 
design feedback compensation for minimum settling time in response to a 
step input, and minimum steady state error for a ramp input. 

Using the program to determine the minimum settling time values of 
Ka, and Ka, as a function of T, a family of curves, with parameter Z, 


U 2 


can be constructed on the Ka, vs. T and Ka, vs. T planes. These families 
of curves are shown by figures (22) and (23). Of these two families of 


curves, the minimum settling time Ka, vs. T plane is by far the most 


I 
interesting. The minimum settling time Ka, vs. T loci was shown by fig- 
ure (19) to be hyperbolic. Furthermore, section 3.3 outlined the presence 
of discontinuities on parameter planes for systems with more than two 
parameters. These discontinuities again appear on the minimum settling time 
parameter planes, figures (22) and (23), and are best demonstrated by fig- 
ure (23). Inspection of figure (22), the minimum settling time Ka, vs. T 
plane, also reveals these discontinuities on the loci of minimum settling 
time for the cases of Z less than zero (right-half s-plane zeros). 

It has been demonstrated, [2] and [7], that for a type one sampled- 
data control system, the steady state error in response to a ramp input is 
inversely proportional to the plant root locus gain; similiar to the "velo- 


city lag'’ error in continuous systems. Consequently, to minimize the ramp 


error and at the same time design for minimum settling time response to a 
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step input, the designer must maximize the parameter K (or Ka, when using 
unity position feedback) on the specified line of constant sampling in- 
terval T = 0.5. Inspection of figure (22) shows the importance played in 
this selection of the plant gain by the discontinuity phenomena. From 
figure (22), the maximum Ka, attainable for T = 0.5 is Ka, = 12.7, when 

Z =-1.5. For Z less than -1.5, the discontinuity asymptote lies to the 
left of T = 0.5, and a negative value of Ka, (positive feedback) is re- 
quired for minimum step input settling time. For Z greater than -1.5, 
lower values of Ka, are required for minimum settling time for a step in- 


put and therefore a larger ramp error is experienced. 


Therefore, the design criteria of this example is achieved for: 


K = 12.7 
ay = 1.0 
ay = 0.656 

Z= -1.5 
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CHAPTER V 
CONCLUSIONS 

The application of Mitrovic's Method to feedback compensation of 
sampled-data control systems has been shown to be quite useful in the 
minimization of system settling time. The parameter plane for the second 
order system presents an easily obtainable picture of the system stability 
and settling time as a function of the system parameters. Using the 
methods described herein, the second order parameter plane development 
results in a two parameter analysis method requiring about the same total 
effort to construct as conventional one parameter methods. 

In the case of higher order systems, the effectiveness of the Mitrovic 
Method as a design tool is somewhat reduced. Foremost, the plotting of the 
various loci requires at least a desk calculator, preferably a digital 
computer. If the number of system parameters can be reduced to only two, 
then the parameter plane is still a good tool for minimizing settling time. 
Unfortunately, the determination of the minimum settling time compensation 
from the higher order parameter plane requires a more tedious study of the 
plane than for the second order case. The intersections of all the loci 
of constant Gz and constant GZ must be examined in detail in search of 
the circle of minimum radius in the z-plane encircling all the closed loop 
roots. 

For higher order systems with more than two degrees of freedom, 
analysis with the digital computer is most effective, indeed necessary. 

The vector-matrix state space equations representing the given system 
were shown to lead directly to the determination of the closed loop root 


positions as a function of up to n parameters, n being the order of the 
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system. The strength of this method is the ability to reduce the design 
criteria to a function of almost any physical parameter or condition. 
However, since a large fraction of the feedback compensation problems are 
not more than two parameter in nature, the Mitrovic Method can be an 


exceedingly important tool to the design engineer. 
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APPENDIX I 
STABILITY ANALYSIS OF VECTOR-MATRIX DIFFERENCE EQUATIONS 
The general form of the state space vector-matrix equation describ- 
ing a linear system with a sampled input is, 

X{(k + 1)T] = FX (kT) = A u(kT). (I-1) 
In this form, the states of the system, K, are defined as the control 
system output and its first n-1] derivatives, where n is the order of the 
linear system. For a feedback control system, the control vector, u(t), 
is the difference between the command input, R(t), and the sum of the 
feedback quantities, ie; 

u(t) = R(t) - AY x(t), (1-2) 
where he is the transpose of A, the coefficient vector controlling the 
magnitude of feedback for each state variable. In discrete time form, 
equation (1-2) is written as 

u(kT) = R(kT) - A X(kT) . (1-3) 
By substituting equation (1-3) into equation (1-1), the familiar vector- 
matrix difference equation for a sampled data control system is obtained 
as 

K{(k + 1)T) = FXCRT)-Aa® xcer) + AR(KT) , eS 
or by imbeding the feedback quantities directly into the F matrix, 
X((k + 1)T) = & X(kT) + A R(KT) (1-5) 
With the feedback thus imbeded in the transition matrix, - , the 
stability analysis of a given system can be confined to an analysis of the 
partial state equation, 
X[(k + 1)T] = X(kT) Gee) 


which is an initial value problem. If the state trajectory returns to a 
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given equilibrium, or singular point, in the state space for any initial 
state vector, it will also move to a predictable translation of this singu- 
lar point and achieve equilibrium for the input R(t) which can be described 
by the initial state vector. The dynamic stability of the system is thus 
ensured. 

Hahn, [1], studied the stability of equation (1-6) using Liapunov's 
Direct Method, and concluded that the stability of equation (Ie6) was 
ensured when the logarithm of the eigenvalues of the ¢ matrix were all 
less than zero in magnitude, that is 

Ind; < 0 (1-7) 
This is equivalent to 

In (Re A, + j Imi) < 0 
Or In ([); | ae )  <aee 
or ha Ag + j ¢ D + 2nTT )< 0 
which is satisfied if and only if | 

ai 
Immediately, a direct analogy of the unit circle in the yD - plane can be 
drawn to the unit circle of the z-plane. This analogy can also be demon- 
strated by considering the z-transformation of equation (I-6) as: 

2X(z) = & X(z) + 2X(0) 
(2I - S )X(z) = 2X(0) 

Or, X(z) = 2(2I - £ ae X(0) 
Now, the inverse of the matrix (zI -@) is the transpose of the cofactor 
matrix, divided by the determinant value of (2I - & y- or, if the cofacuer 


transpose matrix is denoted as [C] 
-] [Cc] 
a - SS eR en 
Me f ) {zl - al 
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Since 2 is just a scalar variable over a complex region, the determinant 





X(0) 


{ 21 - G | is nothing more than the characteristic equation of the 
matrix, whose characteristic values or eigenvalues are the closed loop 
root locations in the z-plane, or ,A -plane. Similiarly, in equation 
(I-11), the eigenvalues of the F matrix are the open loop root locations 
in the z-plane ( r -plane). The vector-matrix formulation is ideally 
suited for digital computation of the state variable time response in 
discrete form. Even more important, it is also suited for determination 
of the system dynamic stability by digitally locating the eigenvalues of 


the $ matrix as a function of any variable in the $ matrix. 
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APPENDIX II 
FORTRAN PROGRAM FOR MINIMUM SETTLING TIME PARAMETER PLANE 


P2OGRAM DESIGN 

PROGRAMMER Le Die NACE 

THIS PROGRAM COMPUTES AND GRAPHS THE VALUES OF THE 
PARAMETERS KA1 AND KA2 REQUIRED FOR MINIMUM PROTOTYPE 
RESPONSE OF ANY SECOND ORDER SAMPLED-DATA FEEDBACK 
CONTROL SYSTEM AS A FUNCTION OF THE SAMPLING INTERVAL Te 
INPUT DATA IS . 
CARD 1 VALUES OF Ps» THE NEGATIVE OF THE PLANT S=PLANE 
POLE LOCATION IN DECREASING MAGNITUDE ORDERs FORMAT 2F10e5 
CARD 2 COEFFICIENTS OF THE S=PLANE NUMERATOR OF PLANT 

IN ASCENDING ORDERs FORMAT 3F10-25 

CARD 3. T MINe T MAXs DELTA Ts IN FORMAT 3F100e5¢ T MIN IS 
THE MINIMUM VALUE OF T SCANNEDe T MAX IS THE MAXIMUM 
VALUEs AND DELTA T IS THE INCREMENTAL VALUE OF T USED IN 
CALCULATIONS. MAXIMUM NUMBER OF POINTS COMPUTED IS 500. 

T MIN MUST BE GREATER THAN ZERO. | 

CARD G4 FIRST LINE OF THE GRAPH TITLEs COLUMNS 1-486 

CARD 5 SECOND LINE OF THE GRAPH TITLEs COLUMNS 1-48 

CARD 6 EXSCALEsYSCALE FOR GRAPH OUTPUT IN FORMAT 2F10e5 
DIMENSION P(2)98(3)sANS1(500) sANS2 (500) 9 TT(500) 

DIMENSION ITITLE(12) he 

READ 1009P 3 ‘. 

FORMAT (2F10¢5) 

FORMAT (/////) 

FORMAT (//) 

PRINT 75 

PRINT 101 

FORMAT (47HVALUES OF P» NEGATIVE OF S=PLANE POLE LOCATIONS) 
PRINT 76 3 

PRINT 100,5(P(I)sI=192) . 

READ 10298 : a 
FORMAT (3F10e5) 

PRINT 75 

PRINT 103 

FORMAT (47HPLANT NUMERATOR COEFFICIENTS IN ASCENDING ORDER) 
PRINT 76 

PRINT 1025(8(1) sI=193) 

READ 102s TMINs TMAXsDELTAT 

PRINT 104 : 

FORMAT (/////53X95HT MINs 6X 95HT MAX ¥6X97HDELTA To//) 
PRINT 1025 TMINsTMAX sDELTAT 

FORMAT (6A8) 

READ 2006 (ITINUEt1 II =e) 

READ 200s(ITITLE(I)s1=791 12) 

PRINT 201 

FORMAT (/////912HGRAPH TITLESs//) 
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PRINT 200s(ITITLE( I) 51=7912) a 
READ 2025 EXSCALEsYSCALF a 

202 FORMAT (2F1005) te 
PRINT: 203 in 
202 FORMAT (/////910X%s THX~SCALE 9 10X 9 7THY=SCALE 9//) Vi 
PRINT 204 ,EXSCALE»sYSCALE ih 
204 FORMAT (8X9F100599XsF1065) Hu 
LAST = 0 : (Vi) 
B= P (9) r : We 
Ee =P ian ss te 
BO= 5 te i 
B1=B(2) (y\ 
B2=6(3) ih 
0 a 
T=TMIN 7 
195 K=K+1 i 
fe) (eden meri. 1ie I) 
116 9U=1./EXPF(P1*T) nt 
V=1le/EXPF(P2*T) Mi 
PH11=P2*U/(P2-P1)+P1*V/(P1-P2) | aN 
PH12=U/(P2-P1)4+V/(P1-P2) ey 


PH?) =(USV=1. Ie CUP 1eP2) /(PIFe2 Noe 
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HOO 1G (N=V) estes eZ ; yt 
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1+((B2*P2*P2—-6 142780) 7 (P24P2—P Eee Zoek 
DO2=((681=P1=P2)* (Pi) 80-4 eA Ae eee 
f+((B1=-P1=P2)*(=P2y 20 -P ley 7 ee 
Go 10, Tie 


Se Sak 
SOF SY es. 
a f 


oe 
—— 







ANS2(K)=DUMA2/DENOM 

(Ck aa 

IF (TMAX—T) 11551159114 
114 T=T+DELTAT 


111 D1l=((B1+B0*T)*P1-BO) /(P1¥*P1)+((B2*P1*P1— Ap 1480)%U1/1 Pee NV 
D2=B0O/P1~—-((B2*P1-B1)*P14+R0) *U/P 1] te 
ea dae Hi 

112 D1=82+B1*T+BO*T#T/26 : NG 
D2=B1+B0*T Wis 
PH1l11=l. it 
PH12=T fo) 
PH21=0. , 
PH22=1. iW 

113 DENOM=D1* (PH21*D1-PH11*D2)-D2*(PH12*D2=PH22*D1) ne 

 -DUMAT=(PH] 1+PH22 ket 2 1b Pia ae aa is 
1+D2* (PH11*PH22-PH12*PH21) ; qi 
DUMA2=-D1*(PH11*PH22~PH1] 2*PH21)—-(PH11+PH22) *(PH12*D2—-PH22*D1) ti 
ANS1(K)=DUMA1/DENOM nt 
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115 PRINT 116 
L16 FORMAT 


(TTC I) sANSL (J) 9ANS2(0J5) sJ=1 9K) 


L117 FORMAT(5X»F100¢5910X9F1005910X9F 1065) 


PRINT ll1l7s 


(KsTTsANS1 909094HKAI1 


16929% 09859851 sLAST) 


CALL DRAW 


CALL DRAW, (Ke TT »sANS2909094HKA2 


1099909098589] LAST) 


END 
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